Entanglement, purity and energy: Two qubits vs Two modes 
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We study the relationship between the entanglement, mixedness and energy of two-qubit and two- 
mode Gaussian quantum states. We parametrize the set of allowed states of these two fundamentally 
different physical systems using measures of entanglement, mixedness and energy that allow us to 
compare and contrast the two systems using a phase diagram. This phase diagram enables one 
to clearly identify not only the physically allowed states, but the set of states connected under an 
arbitrary quantum operation. We pay particular attention to the maximally entangled mixed states 
(MEMS) of each system. Following this we investigate how efficiently one may transfer entanglement 
from two-mode to two-qubit states. 

PACS numbers: 03.67.-a 



I. INTRODUCTION 

In this paper we present a parametrization of the states 
of two qubits and the Gaussian states of a two- mode con- 
tinuous variable system in terms of entanglement, purity 
and energy. This allows one to compare these two funda- 
mentally different physical systems using a entanglement- 
purity-energy (EPE) phase diagram. The relationship 
between these three properties gives us some insight into 
the set of allowed states in both systems and, in par- 
ticular, the distinction between separable and entangled 
states, both pure and mixed. Also there has been consid- 
erable interest in recent years in the maximally entangled 
mixed states (MEMS) of two qubits and two modes. Here 
we re-examine this type of state and its dependence on 
the energy in the system. 

Entanglement is the resource used in many quantum 
information protocols including quantum teleportation, 
quantum cryptography and quantum communication. In 
the simplest form, these protocols ideally require maxi- 
mally entangled pure states. It is highly likely though 
that in any real implementation the states will be mixed 
to some degree. It is then desirable to know how this 
mixedness limits the amount of entanglement. 

Qubits have been the mostly studied basic informa- 
tion unit in quantum information theory. As a two-level 
quantum system they have been experimentally realised 
in many different physical set-ups. Two qubits provide 
the simplest discrete quantum system where non-trivial 
properties of entanglement can be studied. The entan- 
glement can be unambiguously quantified using different 
measures and the maximally entangled pure states are 
the so-called Bell pairs. For mixed states, those states 
which maximise the entanglement for a given degree of 
purity are the MEMS states introduced by Ishizaka and 
Hiroshima 1] . They were studied in more depth by Ver- 
straete et al. in Ref. 0] who derived the unitary oper- 
ation which must be applied in order to maximise the 



entanglement of a state p. An Ansatz for the MEMS, 
which is seen to be locally equivalent to a diagonal state 
transformed under one of the aforementioned unitary op- 
erations, was presented in Ref. ; 3j. Different measures 
of entanglement do not have the same dependence on a 
given measure of mixedness. Thus the MEMS depend on 
which measures are chosen and this fact is investigated 
in detail in Ref. 

The previous studies of maximally entangled mixed 
states use entanglement and mixedness to characterize 
the set of two-qubit states. Indeed, in Ref. ,4j the au- 
thors use a concurrence-purity (CP) phase diagram to 
graphically present the set of physically allowed states. 
A physical realization of two qubits could be two two- 
level atoms or two polarized photons where the | 0) and 
1 1) qubit levels are the ground and excited atomic lev- 
els in the former or vertical and horizontal polarization 
states in the latter case. It is the case that the energy 
in the system has a direct influence on the entanglement 
and mixedness properties. For instance, if we know there 
are zero or two excitations in the system composed of 
two atoms we can tell without ambiguity that the two- 
qubit is both pure and separable. For other energies the 
connection between entanglement, mixedness and energy 
is not so clear. The first aim of this paper is to explore 
this relationship with particular emphasis on "frontier" 
states - those states that lie on the edge of the physically 
allowed states. 

We will also investigate the relationship between en- 
tanglement, mixedness and energy for two- mode Gaus- 
sian states. Interest in the information processing abil- 
ities of continuous variable quantum systems has grown 
dramatically in recent years due mainly to the experi- 
mental benefits of generating and measuring entangled 
states compared to a discrete quantum system. The sim- 
plest system one can then study is that of two modes 
of, say, an electromagnetic field. Also, by restricting to 
Gaussian states, those states with a Gaussian Wigner 
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phase space representation, there exists solid measures 
of entanglement when the energy in the system is finite. 
The idea of a Bell state in this system is not clearly de- 
fined as it is possible in theory to have arbitrarily large 
entangled states limited only by the ability to squeeze 
states. 

Extremally entangled states of two-mode Gaussian 
states were studied in Refs. 0, |^ where it was demon- 
strated that by fixing not only the global purity, but 
also the purities of the reduced density matrices of both 
modes, there exist maximally and minimally entangled 
states. In what follows we will simply fix the total en- 
ergy in the system and then look for the physically al- 
lowed states at this energy, including those that are max- 
imally entangled. The EPE phase diagram offers a nice 
visual aide to the allowed states. The distinction between 
separable and entangled states is clear to see and, while 
maximally entangled states coincide with the previous 
work, in this picture the minimally entangled states are 
just the separable states. 

The rest of the paper is organised as follows. In section 
[n]we cover the case of two qubits in some detail with par- 
ticular regard to the extremal states. Section ITU covers 
the two-mode Gaussian states where we find some sim- 
ilarities and differences with previous efforts to describe 
these states. In section llVl we look at how entanglement 
can be transferred from two-mode to two-qubit states. 
We conclude in section with some closing remarks. 

II. TWO-QUBIT STATES 

Two qubit states are completely specified by 15 pa- 
rameters and as such, it is difficult to find an appropriate 
parametrization in which to view the state space. There- 
fore we want to reduce the number of parameters involved 
while simultaneously retaining as much non-trivial infor- 
mation about the states as possible. 

We choose to use the entanglement, purity and energy. 
First we need to define the quantities we will use to quan- 
tify the entanglement, mixedness and energy for qubits. 
For entanglement we choose to employ the concurrence. 
This is an easily calculable measure for the entanglement 
of two qubits in the state p given by the explicit formula, 

CO) = 2 max{/ij} - ^ to > W 

i 

where {pj} are the square roots of the eigenvalues of the 
matrix R = p(o y ® a y )p*(a y ® a v ) and p* is the complex 
conjugate of the density matrix p. C{p) ranges from zero 
for separable states to one for maximally entangled Bell 
states. The tangle, 

r[p) = [C(p)} 2 , (2) 
and the entanglement of formation, 



where x± = (1 ± y/l - [C(p)] 2 )/2 and, 

h(x) — — xlog 2 x (4) 

and the concurrence are equivalent measures of entangle- 
ment. The negativity, 

llp r ll - f 

N{p) = ^ , (5) 

and the logarithmic negativity, 

Ljv(p) = log 2 ll/H , (6) 

where = TrvAAt and p r is the partial transpose of 
p with respect to one or other of the subsystems, are two 
measures of entanglement we could also have used here. 
The ordering of states is changed depending on the choice 
of entanglement measures jallQ. In Ref. the au- 
thors investigate in detail how the maximally entangled 
mixed states depend on the choice of measures for both 
entanglement and mixedness. For instance, when purity 
is chosen to measure the mixedness and concurrence the 
entanglement, the MEMS are those states defined in Eq. 
(|TS|l whereas if the negativity is used to measure entan- 
glement, the Werner states, 

p w = r\^\ + ^^-I , (7) 

where \ ^>) is a maximally entangled state and I is the 
identity, achieve the highest degree of entanglement for 
a given amount of purity. Here we are interested in what 
parameters we can use to view the set of physically al- 
lowed states of a two-qubit system in such a way as to 
gain as much information as possible. We take gener- 
ally the entanglement, the mixedness and the energy, 
while the measures we choose to quantify these prop- 
erties are not so important. It only affects the shape 
of the entanglement-energy-purity diagram as can seen 
when comparing Fig. 3 to Fig. 4(a). 

To quantify the mixedness of a state p we will use the 
purity, 

Pip) = Tr[p 2 ] • (8) 

In general, for d-dimensional systems Pip) ranges from 
1/d for completely mixed states to 1 for pure states. It 
is closely related to the linear entropy measure of mixed- 
ness, S'l(p) = jzt(1 — P(fl))> d > 1, which ranges from 
zero for pure states to one for the maximally mixed d- 
dimensional state. 

The energy of two two-level atoms with resonance fre- 
quency ujq is defined here as the expectation value of the 

Hamiltonian, H atom = Thoqcj'P + hw^a^ , where a z is 
the Pauli matrix satisfying <r z \j) = (— l)^ 1 ^! j), for the 
ground and excited atomic states {| 0), 1 1)}. We set h 

and too to 1 leaving H = + ai 2 \ We will use the 
average excitation number in the system, 



EoFip) = h(x + ) + hix_) , (3) 



E{p) = 1 + Tr[Hp] 



(9) 
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E(p) ranges from zero (for p = | 00) (00 | ) to two (p = 
| 11) (11 1) for the physically allowed two-qubit states. 

A concurrence and purity phase diagram is one method 
which tells us much about the set of physically allowed 
states. This phase diagram is shown in Fig. 1. It is 
bound by the maximally entangled mixed states defined 
later in Eq. i|18|) . Also shown in this figure are the 
Werner states from Eq. (7J. We now add another param- 
eter the energy of the states. The relationship between 
the energy, purity and entanglement is already known 
for certain extreme two-qubit states. These are the pure 
states, the separable states and the maximally entan- 
gled mixed states mentioned above. The entanglement- 
purity-energy (EPE) phase diagram for these states is 
shown in Fig. 2. Each point in these diagrams repre- 
sents the set of states, 



Thus we first consider the following density matrices, 



{p:E(p) = E,C(p)=C,P(p) = P} 



We first go through explicitly the boundaries of Fig 
which are more clearly shown in Fig. 4(a). 



Pure States 



(10) 
2 



f Pi 

p 2 

p 3 

pi 



(13) 



with E(p) = 1, C{p) = 0, P(p) = 2p\ +p 2 2 +pl and 2 Pl + 
Vi + P3 = 1- States of this form include the completely 
mixed state but the pure states | 00) or | 11) cannot be 
expressed in this way. We shift the energy of the state p 
by a fixed amount, e, to get p e , 



(14) 



so that the energy is now E(p e ) — 1 + e and the purity is 
P(pe) = (pi+ e ) 2 + (P2~ e ) 2 +P3+Pi- The minimum purity 
for a given energy, E, is thus P m in = z(l + 2(-E — l) 2 ) for 
^ < E < | to preserve positivity of the density matrices. 
In order to calculate the minimum purity of states p with 
energy E(p) > | we take the following density matrices, 




For the pure states, P — 1, we have, 

\*)=a\a 1 )\0 1 )+b\a 2 )\fa) , 



(11) 



by the Schmidt decomposition, where a, b £ R, < a, b < 
1 with a 2 + b 2 = 1, and (ay | a k ) = S jk = {f3j \ (3 k ). The 
concurrence for | *) is C(|*)(*|) = 2ab. The maxi- 
mum/minimum energy of | ty) for a given concurrence is 
easily seen to be E max / min = 2b 2 . Thus for a concurrence 
C, the pure states on the EPE diagram are contained in 
a circle centered on E = 1,C = 0, P = 1 with radius 1, 
(E— l) 2 + C 2 < 1. This is the blue curve in Fig. 4 (curve 
satisfying P = 1). 



B. Separable States 

Separable states of a bipartite system may be expressed 
as a convex sum of product states, 



Psep 



pf 



(12) 



P = 



( e 

p- e 

1 -p 

\ 



(15) 



The minimum purity in this case is simply, 
P min (p) = ^E 2 ~2E+1 , 



where p — | 00) (00 | when e = p = 1. For these states we 
haveE(p) = 1+eand purity P{p) = (p-e) 2 + (l-p) 2 +e 2 . 

(16) 
(17) 

These three sections are symmetric about E = 1 as one 
would expect and P m in — § f° r each of the relevant ex- 
pressions at E = | and E = |. This is the green curve 
in Fig. 4 (curve satisfying C = 0). 



with < E < \. Similarly, for § < E < 2, 



(p) = -E 2 -AE- 



This form is not particularly useful as it is difficult in 
general to find such a decomposition. Here we want to 
investigate the dependence of the purity and energy for 
the separable states. In particular we would like to know 
the bounds on these states i.e. given a separable state 
p with purity, P, what is the maximum energy p can 
have? This is equivalent to finding the minimum purity 
of separable state with a certain purity. Given the Hamil- 
tonian we use to calculate the energy, it is reasonable to 
assume the separable states with maximum or minimum 
energy for a given purity will be diagonal in the basis 
{| 00), | 01), 1 10), 1 11)}. ' 



C. MEMS 

The maximally entangled mixed states (MEMS) with 
respect to concurrence and purity measures are conve- 
niently expressed using their concurrence, C, as a pa- 
rameter. 

PMEMS = C|0+W+| + i|Ol)<Ol| 

I_£)(|00)(00| + |11)(11|) , 
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for C G [0, §] while, 

Pmems = C\ (f> + )((f) + | + (1 — C)\ 01)(01 1 , (18) 

for C G [|,0], where \<j> + ) = ^(|00) + 1 11)). As they 
are written above the MEMS have zero energy while the 
dependence of purity on the concurrence is given by, 



1 C 2 
P =3 + T ' 



for C G [0, §] and, 



P = C A 



(1-Cf 



(19) 



(20) 



for C £ [|, 1]. We can find the maximum energy these 
states attain by applying local unitary operations, U G 
SU(2)®SU(2). Take U = !7i(0i,ni) <8> f/ 2 (02,n 2 ), with, 



Uj = exp [—i6j{fij ■ a) 



(21) 



'i — ™v l '"JV'J 

where |nj| = 1 for j = 1,2 and <? = (c x , CTj,, ct z ) are the 
Pauli matrices. Applying U gives p' = U pmemsU' with 
energy, 



E(p') = 1 + 



1 



sin 2 2 



(1 - nl,J - sin 2 0i(l-n 2 )] 



(22) 

Thus when C < § the MEMS have an energy in the range 
and | < E < §. For C > |, the same calculation gives 
an energy for the locally transformed states, 

E(p') = 1+(1-C) [sin 2 2 (1 - nl z ) - sin 2 X (1 - n? )] , 

(23) 

and so C < £7 < 2 - C for these states. The MEMS lines 
with maximum and minimum energy are the red curves 
in Fig. 4 (curves beginning at E = 1, P = 1, C = 1). 

The entanglement-purity-energy phase diagram is 
plotted in Fig. 2 for 100,000 randomly generated den- 
sity matrices confirming the boundaries described above. 
There are some points on this diagram worth highlight- 
ing, 

• {C, P, E} = {1, 1, 1} represents all the Bell states, 
1 



±\ 



V2 
1 

71 



(1 00) ± 1 11)) , 
(I 01) ± 1 10)) , 



(24) 



• {1,0,1} represents all pure states of the form, 

\ip)=a\ 00) + b\ 01) + c| 10) + a\ 11) . (25) 

• {1,0,0} represents | 00). 

• {1,0,2} represents | 11). 

• {0, j, 0} represents the completely mixed state, p = 

• {0, \ < P < |, 1 - < E < 1 + ^} represents 
all the separable states, including the completely 
mixed state, which cannot be entangled by unitary 
operations. 



III. TWO-MODE GAUSSIAN STATES 

We now consider a two-mode continuous variable sys- 
tem occupying a Hilbert space which is the tensor prod- 
uct of two Fock spaces. The creation and annihilation 
operators for the jth mode will be denoted by a] and a,j 
respectively and these give rise to the quadrature phase 
operators Xj = (a,j + aj) and Pj = (aj — cfy/i. We will 
only consider Gaussian states of this system which are, 
by definition, states which have a Gaussian characteristic 
or quasi-probability function in phase space. For these 
states it is only necessary to specify the first and second 
statistical moments of the quadrature operators. Denot- 
ing X = (xi,pi, X2,P2), to specify p we need the mean 
values ((A)) and the covariance matrix a where, 



1 



{{X i X j ) + {X j X i ))-{X i ){X j ) 



(26) 



Unitary operations on p correspond to symplectic oper- 
ations on a [3. Thus, 

UpU^ & SaS T . (27) 

where S G Sp(2, (R)) or SflS T = fi and, 



uj 
uj 



1 

-1 



(28) 



A local unitary operation Ui oc = U\ % J7 2 correspond to 
the local symplectic operation Si oc = Si © S2 and each 
Sj satisfies SjiuSj . 

We can set the first moments of X to zero by means 
of local unitary operations which do not affect the en- 
tanglement or purity properties of p. It has been shown 
that by local unitary operations any CM a of a two-mode 
state can also be brought to the so-called standard form 

nana, 



= S t <tS = 



a c + 

a c_ 

c+ b 

c_ b 



(29) 



Writing <j s f = a (we drop the sf subscript from now 
on for clarity) as 



a sf = S T aS = ( 7 



(30) 



the coefficients a, b, c + and c_ arc determined by the local 
symplectic invariants, 

Det a = a 2 , Det = b 2 and Det 7 = c + c_ . (31) 

The global symplectic invariants, invariant under Spu^ 
operations, are, 

Det a — (ab — ci_)(ab — c_) , 

A = Det a + Det (3 + 2Det 7 , (32) 
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the latter being known in the literature as the seralian 
12] while the former is related to the purity by, 



P = 



1 



VDet a 



(33) 



These local and global invariants determine all the entan- 
glement and mixedness properties of a two-mode Gaus- 
sian state p. 

Positivity of p and the commutation relations for 
quadrature phase operators mean that the CM a is re- 
quired to satisfy [Tj, [14J , 



cr + iCl > 



which in turn means, 



A < 1 + Det a , 



(34) 



(35) 



in order to be a proper covariance matrix. Any covariance 
matrix can be written as |15|. 

a = S T vS , (36) 

where S £ £p(4,ij) and, 

v = &\&g{v-,v-,v + ,v + ) , (37) 

is the covariance matrix corresponding to a tensor prod- 
uct of two thermal states. The state v in Eq. I|37[) has a 
density matrix which is the tensor product of two thermal 
states. 



with, 



o {±) - 

Pth — 



P = P { th ® p[t } > 



2 ^/V ± _lN' 



(38) 



liXil , (39) 



and the average photon number in each of the thermal 
n (±) = ^ ■ (40) 



states, p[h' , is 



The quantities f_} are called the symplectic eigen- 
values as they are invariant under symplectic operations. 
In terms of the global invariants Deter and A the sym- 
plectic eigenvalues have the simple expression, 



2z/| = A ± VA 2 - 4Det<7 
Taking v_ <v +1 Eq. (j3"4*|) reads, 
v- > 1 . 



(41) 



(42) 



We can now carry out a similar investigation of the 
relation between the entanglement, purity and energy 
of two-mode quantum states as we have shown for two 
qubits in the previous section. This is the natural way to 



approach the concept of maximally entangled states in 
continuous variable systems because, as is well-known, 
any measure of entanglement can have arbitrarily large 
values for a given purity. By restricting to a fixed energy 
we can compare continuous and discrete level systems 
on the same terms. First we define the quantities under 
consideration in the case of two-mode Gaussian states. 

The quantification of entanglement of Gaussian bipar- 
tite states has been studied extensively in recent years. 
One way has been to utilize the partially transposed 
state (PPT) or Peres-Horodecki condition for continuous 
variable systems [ljj. The PPT criterion is necessary 
and sufficient for the separability of two-mode Gaussian 
states and is defined as the transposition operation per- 
formed on only one of the subsystems. In a two- mode 
CV system this amounts to a mirror reflection in one of 
the momenta operators, 



P -* P=>Pi 



-Pi 



(43) 



for some i = 1, 2. In terms of local the symplectic invari- 
ants this means that, 



Det 7 — > —Det 7 



(44) 



And as this operation does not preserve commutation 
relations, hence it is not symplectic, the global invariants 
A and Deter change to, 

A = Det a + Det ~ 2Det 7 
Deter = Deter , (45) 

for the new CM a corresponding to the transposed den- 
sity matrix p. The PPT criterion thus means that the 
smallest symplectic eigenvalue of cr, 



A - VA 2 - 4Detcr 



must satisfy, 



£>_ > 1 



(46) 



(47) 



In fact, this eigenvalue completely quantifies the entan- 
glement through its relation with the negativity, 



Af = max 



0. 



1 



2£_ 



(48) 



introduced in Ref. and applied to CV systems in 
Ref. 0, and the logarithmic negativity, 



&n{p) — niax{0, — ln£>_} 



(49) 



which arc both decreasing functions of v- , thus quanti- 
fying the violation of Eq. 147[) . 

Recently, Wolf et al. 48j have defined the entangle- 
ment of formation (EoF) for Gaussian states and, fur- 
ther, have shown how to explicitly calculate this quantity 
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for a given CM a. It is an entanglement monotone un- 
der local Gaussian operations and classical communica- 
tion. This quantity has been further studied in Ref. 
and there it is shown that the entanglement of formation 
leads to a different ordering of the states compared to 
that induced by the negativity above. We could use the 
entanglement of formation here but, apart from the sym- 
metric two-mode Gaussian states, there is no analytical 
expression for general Gaussian states. In fact for arbi- 
trary two-mode Gaussian states, finding the EoF is quite 
a non-trivial task. The solution has a ni ce g eometric de- 
scription which is discussed in Refs. [l8l.ll9j ] Ultimately, 
though, it is determined via numerical minimisation so 
we choose to deal with the logarithmic negativity due to 
its more analytical nature. 

To characterize the mixedness we use the purity again 
which for two-mode pure states is simply given by, 



P(a) 



1 



VDct a 



(50) 



For the energy we will use the sum of the average num- 
ber of photons present in each mode, 



E(p) 



(a\ai) 



(4«2) 



(51) 



In terms of the CM, cr, corresponding to the state p, 

E(p) = ^Tra-l + J2(X j ) 2 ■ (52) 

3 

As we have mentioned, the mean values, (Xj), can be set 
to zero by local operations which do not affect the en- 
tanglement or purity. These operations obviously change 
the energy though, so while we set them to zero in the 
following, the energies calculated are the minimum val- 
ues for a given a in the standard form. As an example all 
two-mode coherent states are represented by the vacuum, 
with zero energy. 

The local invariants a and & of a CM a are related to 
these expectation values of the number operators of the 
two modes by 



a = 2ni + 1 , b = 2h 2 + 1 



(53) 



Also a and b are the inverses of the purities of the reduced 
density matrices of the two modes, a = I/Ml = 1/Trpi, 
b = l//i2 = l/TrpI with p\i = Ti^ip. This fact has been 
used to characterize the maximally and minimally entan- 
gled two-mode mixed Gaussian states, for fixed global 
and local purities HJ, where the entanglement is quanti- 
fied using the negativity. 

We will now outline the EPE diagram for the two- 
mode Gaussian states by first examining the extremal 
states. These states have been studied in depth in Refs. 
0, 0| . In Ref. the authors study the extremal en- 
tanglement and mixedness of two-mode Gaussian states. 
They introduce generalized maximally and minimally en- 
tangled mixed states, denoted GMEMS and GLEMS re- 
spectively, which are maximally and minimally entangled 



states for fixed global and local purities. Given how the 
local purities are defined, this would appear to be equiv- 
alent to finding the maximally and minimally entangled 
states for a fixed purity and average number of photons 
in each mode. However, we will see that the minimally 
entangled mixed states for fixed energy are simply the 
separable states. 

In the following we will restrict ourselves to the range 
of energies, < E < 2, where the energy is measured in 
dimensionless units (number of excitations), in order to 
draw parallels with the two-qubit case. 



A. Pure States 

A pure, P = 1, two-mode Gaussian state has a sym- 
metric CM, corresponding to a two-mode squeezed vac- 
uum, with a = b and satisfies c + = — c_ = \J a 2 — 1. The 
minimum energy is simply, 



Emin(p) = a - 1 



(54) 



There is thus a simple form for the symplectic eigenvalue 
of the partially transposed density matrix, 



= a — yo 



1 



(55) 



The dependence of the logarithmic negativity for pure 
states on the energy is shown in Fig. 5 up to a maximum 
energy of 2. Any point on the curve in Fig. 5 represents a 
pure state with CM given by the above values. However 
by applying a local squeezing operator of the form, 



(56) 




we can change the energy of the state to, 
E{U{S sq y P U(S sq )) = acosh2r- 



1 



(57) 



so that the curve shown indicates the maximally entan- 
gled pure two-mode states for the given energy, i.e. the 
two-mode squeezed vacua, and all pure states above the 
curve are physically allowed. 



B. Separable States 

The separable states are all those for which the small- 
est symplectic eigenvalue of the partially transposed CM 
satisfies, 



> 1 



(58) 



or in other words the partially transposed covariance ma- 
trix is still a valid one. The separable two-mode Gaussian 
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states are determined by the value of the energy through 
the condition, 



P > 



1 

ab 



(59) 



which is a direct consequence of keeping a > 0, and im- 
plies that 



P > 



1 



(E + iy 



(60) 



for separable states. The inequality l|60|) comes from con- 
sidering all the combinations of "local" energies in each 
mode which added together give E. The minimum occurs 
when both modes contribute equally i.e. for symmetric 
states. These states are the tensor product of two ther- 
mal states in Eq. i|38|) with equal average photon number 
in both modes. The curve shown in Fig. 6 shows the 
minimum energy a separable two-mode Gaussian state 
can have for a given purity and thus, as above, all the 
physical states lie above this curve. 



GMEMS & GLEMS 



and the purity, 
P 



J{ab-c%){ab-c 2 _) 



(65) 



It is easily seen from Eq. Q46|). and pointed out in Ref. 
[j| , that the smallest symplectic eigenvalue of the CM of 
the partially transposed state, satisfies, 



dv 2 _ 
~dA 



< 



(66) 



E,P 



so that the entanglement of a given state is monotonically 
increasing with the seralian of the CM corresponding to 
the partially transposed state. Thus, in order to find 
the states with maximum entanglement for a given pu- 
rity and energy, we must find where A has a maximum. 
Manipulating Eqs. I|64|) and (|65|l gives, 



A = 4(E+l) 



2\/-^2 +a6(c+ +c_) 2 



(67) 



which is maximised when c+ = — c_ and thus the maxi- 
mally entangled mixed Gaussian states satisfy, 



The Gaussian maximally and minimally entangled 
states for fixed global and localpurities were introduced 
by G. Adesso et al. in Rcfs. 5, 6] through a parametriza- 
tion of two-mode Gaussian states in terms of the local 
purities, 



1 



1 



Mi = 



= 



(61) 



and the global symplectic invariants, the purity and ser- 
alian. By writing the parameters a, b and c± in terms of 
the above parameters, the seralian must satisfy, 



where A min and A max 



are easily calculated. From Eq. 



()46J) it is seen that, 



dv 2 _ 
~dA 



>0 , 



(63) 



so that the entanglement of a state as measured by the 
negativity or logarithmic negativity is a monotonically 
decreasing function of the seralian. Thus the limits im- 
posed by Eq. (|62[) imply that there exist maximally 
and, surprisingly, minimally entangled two-mode Gaus- 
sian states for fixed values of the global and local purities. 

Here we are looking for the maximally entangled mixed 
Gaussian two-mode states for fixed energy. They are 
found by again considering a state, p, with CM a in the 
standard form (|29[) . We now fix the energy of the state, 



E = \{a + b)-l , 



(64) 



A = 4(E+ l) 2 - 



P 



(68) 



These are shown in Fig. 7 for energy 2 (in dimensionless 
units). 

The Gaussian maximally entangled mixed states for 
fixed global and local purities, or GMEMS, as defined 
in Refs. 0, |(| have a simple expression for states with 
a standard form CM. These states saturate the lower 
bound in l|62|) and this means that, 



±\ ab- 



1 

P 



(69) 



(62) and importantly for the definition of entanglement, 



A = 4(E+l) 2 -- 



(70) 



in agreement with Eq. (|68|) . Thus, not unexpectedly, the 
maximally entangled states for a given purity and energy 
coincide with the GMEMS. 

In Fig. 7 one notices the fact that there is no entan- 
glement below a certain purity value. This arises from 
the condition that z>_ > 1 which, together with Eq. 168|l , 
means that, 



P > 



1 



2E+1 



(71) 



in order to be able to unitarily generate entangled states. 
This bound and the condition in Eq. 1|60|) defines a range 
of values of the purity of two-mode Gaussian states, 



1 



(E + l) 



< P < 



1 



2E+ 1 



(72) 



8 



such that, for a given energy E, there can be no entan- 
glement in these states. This is similar to the case of two 
qubits where if the purity is lower than one third, only 
separable states exist independent of the energy. Phys- 
ically, it means that in this range for a given value of 
energy, the states are simply too classically mixed to ac- 
commodate quantum correlations. The range of values 
for the purity increases from zero, at E = 0, reaches a 
maximum when the energy attains the golden ratio, 



E = 



1 + V5 



(73) 



and then falls to zero again as the energy increases to 
infinity. 

It is an interesting coincidence that for the case of E = 
1 both the 2-mode and 2-qubit upper and lower bounds 
agree, \ < P < |. This suggests the maximally mixed 
two-mode Gaussian state, of a given energy E, is that 
state which saturates the lower bound in Eq. H72J1 . Such 
a state would have a CM, 



ed = {E + 1)I 4 



(74) 



where I4 is the 4x4 identity matrix, in analogy with 
the two-qubit case, and corresponds to a state which is a 
tensor product of two thermal states with equal average 
number of photons in each mode as we have seen before 
in Eq. J2§. 

The continuous variable version of the Werner states 
of two qubits in Eq. (0 were introduced in Ref. j2£j and 
have the form, 



Pgw = rpME + (1 - r)pi 



(75) 



where pme is the pure two- mode Gaussian state with en- 
ergy E and pi is the completely mixed state specified by 
the CM in Eq. (|74^l . i.e. the density operator pi is equal 
to a tensor product of two thermal density operators each 
corresponding to a state with the average photon number 
equal to n = E/2 (we remind us that that we consider 
units such that h = uj = 1), i.e., 



Pgw = rpME + (1 - r)p th <8> Pth 



(76) 



where pth is given by Eq. (|39[) in terms of symplectic 
eigenvalues, or in terms of the mean excitation number 
n we have 



Pth 



1 



(n+ 1) ^ 

y ' 3=0 



1 



(77) 



These states are not Gaussian however and thus do not 
appear on our EPE phase diagram. 

The GLEMS are states which by definition saturate the 
upper bound in Eq. 116211 . From Refs. states satisfy 

this condition if and only if the symplectic spectrum is 
V- = l,v+ — 1/P and this leads to, 
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1 

a b 



(A-(a~b) 2 ) 2 -(A-(a + b) 2 Y 



(78) 



for these states described by a standard form CM and 
thus, 



A = 4(E+iy-[l + — 



1 



(79) 



The first term arises out of the fact that the minimum 
value of A occurs when the states are described by a 
symmetric CM. The GLEMS actually describe a range 
of states of which those with the minimum entanglement 
are shown in Fig. 8 for energy 2 and they exist in the 
darker shaded region between the two curves shown. The 
GLEMS and GMEMS coincide when the global purity 
can be written in terms of the two local purities and 
these states are called GMEMMS or Gaussian maximally 
entangled mixed states for fixed marginals. As noted in 
Refs. HQ, the GLEMS contain separable states which 
are readily seen in Fig. 8. The range of purities for which 
these states can tolerate separable states is, for a given 
energy E, 



1 



2E + 1 



< P < 



1 



V2e t +IeTT 



(80) 



We can now put together the pieces of the EPE dia- 
gram for 2-mode Gaussian states. The resulting graph, 
shown in Fig. 9, follows the outline of a quarter of a 
fruit bowl. Now in this picture the idea of minimally en- 
tangled mixed states for a fixed energy has no meaning, 
unless you consider the separable states as such. Thus 
the GLEMS do not arise naturally in our picture. This 
is understandable as we are using only three parameters 
to describe the set of states instead of four and thereby 
lose some information. 

Some points on this diagram worth highlighting are 
the following, 



{C,P,E} = {0,. 



( E+1 yi , E} represents all the com- 
pletely mixed two-mode Gaussian state. 

• {0,1,0} represents all tensor product of the vacuum 
state of two modes 

• {0; (E+iy 1 — P — 2E 1 +i ' ^} represents all the sepa- 
rable states, including the completely mixed state, 
which cannot be entangled by arbitrary unitary op- 
erations. 



IV. ENTANGLEMENT TRANSFER: 2 MODES 
TO 2 QUBITS 

Now we complete our study of the two sets of states, 
those of two qubits compared to the two-mode Gaussian 
states, by considering how efficiently one may transfer 
entanglement and en erg y from a two-mode field to two 
qubits. In Refs. pll I22L I23L |24| entanglement swapping 
between qubits and continuous variables has been investi- 
gated in detail. Here we are particularily interested in the 
role of energy in the entanglement transfer process. We 
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will employ the basic resonant Jaynes-Cummings Hamil- 
tonian |25| . 



(81) 



where aj and a,j are the creation and annihilation opera- 
tors for the jth mode and are the raising and lowering 
operators for the jth atom, to describe the interaction of 
the two modes with the two two-level atoms. The situ- 
ation envisaged is that the first mode interacts with the 
first atom only and similarily for the second mode and 
the second atom with both atom-field couplings equal, 
Ai = A2 = A. We will analyse three different initial 
states of the field while the state of the two atoms will 
always be in the ground state (| g)\ | 5)2) initially. 

First we consider when the field is in a superposition 
of one photon in the first mode and zero photons in the 
second and vice versa, 



1 



|^(t = Q)) = -^(|Q)i|l) a 



1>i|0)j 



(82) 



There has been some recent discussion [26|, |27| related 
to this state. Intuitively this state is an entangled state 
of two modes and should generate entanglement between 
the two qubits. Indeed, as the input state is pure, in this 
case and all the following cases, we can use the von Neu- 
mann entropy of the reduced state of one of the modes 
as a measure of the entanglement of this state, 



S(|V/(0))) = -Tr ft mp,-, i = 1,2 



(83) 



where pj is the reduced density matrix of the first or 
second mode. 

The initial energy in the two-mode field is one exci- 
tation and the entropy is S(\ i/>/(0))) = In 2. When we 
allow the field and atoms interact for a time, t, under the 
Hamiltonian in Eq. (|81|l . the combined output state for 
the two-mode field and the two qubits is, 



|V>/,at(*)> = COs(At) I V/(0)> |<?<7) 

-i sin(At)|0) 1 |0) 2 (l5e) 



\eg)) 



V2 



(84) 



where | 0)j represents the vacuum state of the jth mode. 
Thus at any time t such that gt = ^ the two qubits 
are in a maximally entangled Bell state. At these times 
the initial energy in the field is completely transferred to 
the state of the atoms and the joint state of the atoms 
and two-mode field is separable and pure. What about a 
similar but more general input field state of the form, 



1 



1^/(0)) = -^(|0)i|n) 3 + 



ra)i|0> 2 ? 



(85) 



The initial entanglement is In 2 as previously but there 
are now n excitations present. It turns out that, after 



the same evolution under Hamiltonian in Eq. 1811) . 
resulting two-qubit density matrix is given by, 



the 



Pit) 



f 2 



cos 2 Xt^/n 



V 





sin 2 Xty/n 
t i sin Xty/n 








Xty/n 




S n ,i sin 2 Xty/ii 


°/ 

(86) 

so that for n > 1 no entanglement is transferred from 
the two-mode field to the pair of atoms and their state 
remains separable even though energy is continuously 
swapped between the field and the atoms. 
The next case we consider is the input field, 



I MO)) 



1 



v^Vl + e-H 



= (|0)i|a) a + |a)i|0) 2 ) , (87) 



where | a) is a coherent state. The initial energy in this 
state, N = (Ni + N 2 ), is, 



N 



a 



1 + e-H 2 



(88) 



and the initial entanglement ranges from to In 2 depend- 
ing on the value of a. After evolution under Eq. (|81|l for 
a time t, the state of the two qubits is well-described by 
a density matrix of the form, 



Pit) 




(89) 



with 



-M 



2(l + e"H 2 ) 



sin 2 (At) 



(90) 



The concurrence of this state is C = 2z which is max- 
imised when there is on average one excitation initially 
in the field modes, i.e. N = 1. This is shown in Fig. 
10 where the maximum entanglement generated between 
the two atoms is plotted as a function of the initial energy 
in the field. We have also plotted in Fig. 11 the depen- 
dence of the maximum entanglement generated between 
the two atoms on the initial entanglement present be- 
tween the two modes. It is interesting that for the initial 
state of the two modes in Eq. H87J1 the degree of entangle- 
ment of this state need not be the maximum allowed in 
order to achieve the maximum amount of entanglement 
possible between the two atoms. The fact that the entan- 
glement between the two atoms is maximised when there 
is on average one excitation initially in the field is not so 
surprising given the previous input states considered. In 
any case, the amount of entanglement generated between 
the two atoms is very small. 

Finally we consider the case when the input state of 
the field is a pure two-mode Gaussian state, 



(91) 
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with 7 = tanhr where the above state can be gen- 
erated by applying the two-mode squeezing operator 
S(r) — exp[— r(a\a\ — a%a2)], to the two-mode vacuun 
state. The initial energy in the held is given by, 



27 2 



I-7 2 ' 
and the initial entanglement is, 

Sin — hi(l - 7 2 ) + E in In 7 



(92) 



(93) 



The state of the two qubits after a time evolution is given 
by, 



P(t) 



( a x 

6 

6 

\x d 



(94) 



where, 



6 
d 



^^ 7 2 "(cos2V^At + l) 2 , 

n 

^^ 7 2 "sin 2 2v^Ai , 

n 

iA^ 7 2 "(cos2^At-l) 2 , 

n 

7 2 " +1 (cos 2Vn~+T\t - 1) 

n 

x(cos2 v / nAt + 1) , 



(95) 



and A = (1— 7 2 ). In Fig. 12 the maximum entanglement, 
C ut i generated between the pair of atoms and the purity, 
P, of the state of the two atoms at this time is plotted 
as a function of the initial energy in the field, Ei n . The 
maximum entanglement generated between the atoms is 
quite high, C™" 31 ~ 0.9, in this case but it does not oc- 
cur when there is initially one excitation in the field as 
was the case for the previous input states. In Fig. 13 we 
plot the dependence of the maximum entanglement gen- 
erated on the initial entanglement in the field. Again we 
see that it is not necessary to initially have the most en- 
tangled two-mode state possible in order to generate the 
most maximally entangled two-qubit state possible by 
evolution under the resonant Jaynes-Cummings Hamil- 
tonian. 



V. CONCLUSION 

We have set out to show that the energy of a state is an 
important consideration when specifying the maximally 



entangled mixed states of either two qubits or a two-mode 
continuous variable system. 



For two qubits there is a continuous range of MEMS 
lines on the entanglement-purity-energy (EPE) diagram, 
each depending on the average number of excitations in 
the two-qubit system. These states, together with the 
separable and pure states, form the surface of the volume 
containing all physically allowed quantum states. The 
states in <|18|) are those which maximise the concurrence 
for a given purity and have exactly one excitation. They 
include the Bell states (|24fl . If a state has E ^ 1 then 
it can still be considered a maximally entangled state 
albeit one less entangled than the states (|18f) for the same 
purity. 



For the case of two-mode states it is natural that the 
energy be fixed in order to specify the maxmally entan- 
gled mixed states. Through this parametrization we have 
provided a nice way view the set of two-mode Gaussian 
states. From this EPE diagram the set of physically al- 
lowed states of these two fundamentally different systems 
can be compared. There are some interesting similarities 
we have pointed out which go beyond the fact that both 
systems are completely specified by a 4 x 4 matrix. The 
EPE diagrams are also an aide to view the allowed dy- 
namical evolution of these states. 



Finally we have looked at how entanglement can be 
swapped between these two systems using a simple reso- 
nant Jaynes-Cummings model. We have shown that the 
initial energy and the initial entanglement present in the 
two-mode field are important factors as to how efficiently 
the entanglement is transferred from one system to the 
other. Physically this is a method of distributing entan- 
glement between two distant atomic systems and what is 
interesting is that the two-mode Gaussian states are very 
effective at entangling the two qubits. Given that these 
states are experimentally less demanding to produce, this 
is a nice result. The natural extension of this idea to a 
two-mode field interacting with N atoms may serve as a 
quantum memory [§8| and this will be studied later [2^| . 
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FIG. 1: The concurrence vs purity phase diagram of two qubits. The lower curve, Cw, indicates the Werner states, while the 
upper curve corresponds to the MEMS. All quantities are dimensionless. 
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FIG. 2: (Color online) Numerically determined entanglement-purity-energy (EPE) phase diagram for 100000 randomly chosen 
two-qubit states where the entanglement is measured in terms of the concurrence. All quantities are dimensionless. 
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FIG. 3: (Color online) The boundaries of the entanglement-purity-energy (EPE) diagram for two-qubit states when the log- 
arithmic negativity is used as the measure of entanglement. The MEMS are the Werner states from Eq. (7) in this case 
and form a line rather than a plane as in the case when concurrence is used to measure the entanglement. All quantities are 
dimensionless. 




(d) View of Fig. 4(a) looking 
along the Concurrence axis. 



FIG. 4: (Color online) The entanglement-purity-energy (EPE) phase diagram of the (i) Pure (blue curve) (ii) Separable (green 
curve) and (iii) maximally entangled mixed states (red curves) of two qubits. All quantities are dimensionless. 




FIG. 5: Graph depicting how the energy, E, of pure two-mode Gaussian states depends on the amount entanglement, £ , 
between the two modes. The physical region lies in the shaded region above the curve. All quantities are dimensionless. 




FIG. 6: Graph showing how the energy, E, of the separable two-mode Gaussian states depends on the purity, P, for energy 
values up to 2. The physical region lies in the shaded region above the curve. All quantities are dimensionless. 




FIG. 7: Graph showing how the entanglement, £ , of two-mode Gaussian states depends on the purity, P, for a fixed energy 
value of 2. The upper curve indicates the maximally entangled mixed states for fixed energy and the physical region is the 
shaded area below this curve. All quantities are dimensionless. 
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FIG. 8: Graph showing the GMEMS and GLEMS of two-mode Gaussian states for fixed energy value of 2. The GMEMS are 
those states lying on the upper curve while the GLEMS are those states which lie in the darker shaded region between the two 
curves. All quantities are dimensionless. 



FIG. 9: (Color online) The entanglement-purity-energy (EPE) phase diagram boundaries showing the (i) pure (blue curve 
satisfying P = 1) (ii) separable (green curve satisfying £ = 0) and (iii) maximally entangled mixed (red curves starting on the 
pure curve) Gaussian states of a two-mode continuous variable system for energies, E, in the range < E < 2. All quantities 
are dimensionless. 
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FIG. 10: The dependence of the maximum concurrence, Cut, generated between two two-level atoms by evolution under the 
Jaynes-Cummings Hamiltonian (81), on the initial energy, Ei„, in the two mode field. The field is initially in the entangled 
superposition state of Eq. (87) characterised by the coherent state amplitude a. All quantities are dimensionless. 
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FIG. 11: The dependence of the maximum concurrence, C maX! generated between two two-level atoms by evolution under the 
Jaynes-Cummings Hamiltonian (81), on the initial entanglement in the two mode field. The field is initially in the entangled 
superposition state of eqn. (87) characterised by the coherent state amplitude a. All quantities are dimensionless. 
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FIG. 12: (Color online) The dependence of the maximum concurrence, Cut (lower curve), generated in a pair of two-level 
atoms by evolution under the Jaynes-Cummings Hamiltonian (81), and the purity, P (upper curve), of the state of the two 
atoms at this time on the initial energy, Ei n , of the two-mode field when the field is initially in a pure two-mode Gaussian 
state. All quantities are dimensionless. 




FIG. 13: The dependence of the maximum concurrence, Cout, generated in a two-qubit pair by evolution under the Jaynes- 
Cummings Hamiltonian (81), on the initial entanglement in the two-mode field, Si„, when the field is initially in a pure 
two-mode Gaussian state. All quantities are dimensionless. 



